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NOTE  ON  A NON-TRIVIAL  SIMPLE  EXAMPLE  OF  HIGHER-ORDER 
ONE-DIMENSIONAL  BEAM  THEORY 
by 

E.  Relssner 

ABSTRACT 

One-dimensional  beam  equations  are  derived  for  a rectangular- 
cross  section  beam  consisting  of  shear  webs  and  corner  stringers 
with  one  cross-sectional  axis  of  symmetry.  It  is  shown  that  a 
systematic  statement  of  this  problem  involves  one  more  equation  than 
elementary  theory,  with  this  additional  equation  involving  the  concept 
of  cross-sectional  warping  and  bi-moment.  As  an  application  of  the 
general  results  of  the  paper,  new  conclusions  are  deduced  concerning 
the  location  of  the  center  of  shear  of  the  beam. 
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Note  on  a Non-Trivial  Simple  Example  of  Higher-Order 
One- Dimensional  Beam  Theory1 
by 

E.  Reissner 

Department  of  Applied  Mechanics  and  Engineering  Sciences 
University  of  California,  San  Diego 
La  Jolla,  California  92093 

Introduction.  We  depart  from  the  well-known  fact  that  elementary 
one-dimensional  theory  of  space-curved  beams  involves  a system  of-  six 
equilibrium  equations,  six  strain  displacement  equations  and  six  constit- 
utive equations  for  forces  and  moments,  translational  and  rotational 
displacements,  and  force  and  moment  strains,  and  that  this  theory  may 
turn  out  to  be  inadequate,  for  beams  with  properties  other  than  those 
corresponding  to  compact  and  reasonably  homogeneous  cross  sections.  In 
order  to  gain  insight  into  situations  which  may  be  encountered,  we  consider 
what  we  think  is  the  simplest  non-trivial  example  of  a problem  of  this 
nature,  to  wit  the  problem  of  deriving  the  equations  of  (linear)  one- 
dimensional beam  theory  for  a straight  beam  having  a rectangular  cross 
section,  with  one  axis  of  symmetry,  as  indicated  in  Figure  1.  We  assume 
that  the  bending  stiffness  of  this  beam  is  concentrated  in  four  corner 
stringers  and  that  transverse  shear  and  torsional  stiffness  are  concentrated 
in  four  shear  webs+t.  In  order  to  include  in  our  analysis  the  limiting 
case  of  torsion  of  a thin-walled  open  cross  section  beam,  we  assume  that 
the  four  webs  possess  not  only  inplane  (sheet)  shear  stiffness,  but  also 
transverse  (plate)  twisting  stiffness. 


t A report  on  work  supported  by  the  Office  of  Naval  Research. 

Consideration  of  more  general  configurations  involving  stringers 
and  shear  webs,  but  from  a somewhat  different  point  of  view  may 
be  found  in  |lj. 
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The  relevant  equations  of  statics  of  elementary  beam  theory  with 
which  we  are  concerned  in  what  follows  are  three  equations  for  bending 
moment  M,  transverse  force  Q and  twisting  moment  T,  of  the  form  Q'+q  * 0, 

M'-  Q + m * 0 and  T'  + t * 0.  We  will  show  that  a rational  treatment 
on  the  basis  of  an  elementary  consideration  of  the  three-dimensional 
aspects  of  the  problem,  by  necessity  introduces  a fourth  equation 
of  statics  of  the  form  R'  - P + r * 0.  In  this  equation  the  quantity  R 
has  the  character  of  a one-dimensional  axial  warping  stress  measure  which 
is  effectively  equivalent  to  what  has  been  designated  by  Vlasov  as  "bi- 
moment," and  the  quantity  P has  the  character  of  a one-dimensional  cross 
sectional  shear  stress  measure,  other  than  Q and  T,  the  significance  of  which 
becomes  apparent  in  the  course  of  the  step  from  three-dimensional  formula- 
tion to  one-dimensional  formulation. 

Having  four  one-dimensional  equilibrium  equations  for  Q,  M,  T,  P and 
R we  use  an  appropriate  form  of  the  principle  of  virtual  work  for  the 
derivation  of  strain  displacement  relations  in  a manner  which  is  simple 
but  not  self-evident.  Finally,  we  use  an  appropriate  expression  for  a 
three-dimensional  complementary  energy  density  function  in  order  to  estab- 
lish a system  of  one-dimensional  constitutive  equations. 

As  an  application  of  our  system  of  one-dimensional  beam  equations, 
we  consider  once  more  the  problem  of  determining  the  shear  center  of  the 
cross  section  of  our  beam,  with  this  approach  leading  here  to  more 
detailed  insights  than  are  believed  to  have  been  known  before. 

Statics  and  Virtual  Work.  We  introduce  axial  forces  N. , N0 , shear 
flows  S^,  S2,  Sjj  and  plate  twisting  moments  , M2.  as  indicated 

in  Figure  2,  together  with  load  intensities  n,  s and  m as  also  indicated 
in  the  Figure. 


I 

Denoting  differentiation  with  respect  to  an  axial  coordinate  z by 
primes,  we  then  read  from  Figure  2 the  following  four  equilibrium  equations 

- S2  + S12  + n2  - 0 , Ni  + si  - S12  + nl  - 0 , (la,b) 

2c(S1  - S2)'  + 2c(sl  - s2)  - 0 , (lc) 

2c(S2  + S^)'a  4-  2aS j22c  4-  2c(s2  4-  s^)a  4-  2as^22c 

4-  2[2c(M1  4-  M2)  4-  4aM12l’  4-  2c  (m^  4-  m2)  4-  4am12  - 0 . (Id) 

We  note  that  the  factor  "2"  in  front  of  the  M-terms  in  Eq . (Id)  has  been 
introduced  In  order  to  account  for  the  effect  of  reactive  plate  transverse 
shear  stress  resultants  In  a way  corresponding  to  the  result  obtained 
by  an  explicit  treatment  of  the  problem  of  torsion  in  thin-walled  open 

and  closed  cross  section  beams  (2). 

- 

For  what  follows  it  is  appropriate  to  rewrite  Eqs . (1)  in  terms  of 
one-dimensional  beam  stress  resultants  and  couplee  Q,  M and  T and  in 
terms  of  supplementary  one-dimensional  stress  measures  P and  R,  in  conjunc- 
tion with  analogously  defined  load  intensities  q,  m,  t and  r.  Further- 
more, account  is  taken  of  the  fact  that  the  torque  T is  composed  of  three 

distinct  portions  T^,  T^,  Tq.  The  relevant  defining  relations  are 

2c(S2  - Sx)  - Q . 2c(N1  4-  N2) 

4acS12  " TA  ’ ZcaCSj  4-  S2)  - 

2^Mgds  - 4c (Mj  4-  M2)  4-  8aM12  - 

T “ TA  + TB  + TC  * p * Tc  - TA  • R * 2ca(N2  - Nj)  . (2f,g,h) 

’ 

With  these  we  obtain  directly  from  (lc)  and  (Id)  as  two  of  three  standard 
^equations  of  one-dimensional  beam  theory 


- M , 

(2a, b) 

TC  • 

(2c,<?) 

TB  • 

(2e) 
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Q*  + q - 


0 


T*  + t - 


0 


(3a. b) 


The  third  standard  equation  follows  by  adding  Eq.  (la)  to  Eq . (lb),  in 
the  familiar  form 

M'  - Q + m - 0 

Finally,  the  non-standard  fourth  equation,  involving  the  bi-moment  R and 
the  shear  stress  measure  P follow  by  subtracting  Eq . (lb)  from  Eq . (la), 
in  the  form 

R'  - P + r - 0 , (3d) 

where  use  has  been  made  of  the  defining  relations  (2c),  (2d)  and  (2h). 

We  next  derive  a system  of  strain  displacement  relations  which  is 
consistent  with  the  equilibrium  relations  (3)  by  writing  as  equation 
of  virtual  work 

f (qdv  + u6+  + t«e  + r6*)dz 

- / (Q«>  + Miic  + Ta6ta  + TjjiTg  + Tc«tc  + R6X)dz  . (4) 

Equation  (4)  leads  to  a system  of  virtual  strain  displacement  relations 
upon  eliminating  in  it  the  quantities  q,  m,  t,  r through  use  of  Ecs.  (3), 
and  upon  then  considering  the  resulting  equation  as  an  identity  in  Q,  M, 

Ta.  Tg,  and  R.  In  view  of  the  linearity  of  the  problem  the  step  from 
virtual  to  actual  strain  displacement  relations  is  immediate  and  we  have 
as  the  desired  system  of  relations 


Constitutive  Relations . We  depart  from  the  complementary  energy 


density  expression 


W - 


Nt  NS  cS,  cS~  2aSt~  , 

+ Wi  * + *i®I  + HPU  + 


715- 


(6) 


where  D ■ G t"/6.  We  rewrite  the  various  terms  in  this  expression  in 
s s s 

terms  of  one-dimensional  beam  theory  stress  measures  through  use  of  the 
defining  relations  (2),  as  follows. 


Ni 


N2 

A2E2 


(M-R/a)2  + (W-R/a)  2 
lbc^jEj  16c2A2E2 


M2  . R2  . MR 


(7) 


2aS12 

t12C12 


2*TA 

16c2a2t12G12 


(8) 


cSi 


cS 


2 (Q-Tc/«)2  (Q+Tc/a) 

» — ■ ■ wr'L — m i — rr — r — »«— 


16ct2G2 


ife  |(tpr  *<&)  (,Z  + 7)  + 2(^'e^)t£ 


OT 


7BZ 


_C 

CQ 


(9) 


T2 

- fVfr  * 'Ifa.*-  - 


Dg  " 2 Dgds 
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(10) 


The  meaning  of  the  stiffness  coefficients  D^,  da>  Dq,  Dc 

and  DCq  may  be  read  directly  from  Eqs.  (7)  to  (9).  An  expression  for 
Dp  follows  upon  introducing  the  expression  for  I>s  in  Eq.  (6),  in  the  form 

dB  “ J KGjtJ  + G2t^)c  + 2G12tJ2a]  . (11) 

Having  Eqs.  (6)  to  (10)  we  have  further,  as  a consequence  of 
Castigllano's  theorem, as  the  appropriate  system  of  constitutive  relations 


Y 


TC 

tA 


aw 

" ' 


aw 


a w 


M . R 


X - 


a w 


aw  r . M 

w ^ * 


(12a, b) 
(12c, d) 
(12e,f) 


End-Loaded  Cantilever . Determination  of  Center  of  Shear . We 
assume  as  boundary  conditions  for  the  loaded  end  of  the  beam 

z - 0 ; Q-Qo.T-To,M-R-0,  (13) 

and  as  conditions  for  the  supported  end 


z-L;  v-o,e-o,+-o,f-o.  (14) 

We  solve  the  homogeneous  equilibrium  equations  (3a,b,c),  as  in  the 
elementary  theory 

Q * Qo  * T " T0  » M " z%  • (15a, b.c) 

and  leave  the  homogeneous  equation  (3d)  in  the  form  R'  - 4-  * 0. 

The  constitutive  equations  (12)  in  conjunction  with  the  strain 
displacement  relations  (5)  are  now 


-6- 


t 


(16a, b) 


(16c , d) 

(I6e , f) 


where,  according  to  (2f)  and  (15b)  , TA  + T„  + T,,  « T . 

A D t O 

We  next  introduce  R,  T and  T from  (16)  into  the  remaining  equi- 

v A 

librium  equation  and  obtain  as  one  of  two  differential  equations  for 
t and  0 , 


V - (Da  + Dc)*  + (Da  - Dc)0' 


(17a) 


A second  differential  equation  follows  from  (16)  in  conjunction  with  the 
condition  of  constant  T in  the  form 


<da  + °B  + Dc>«’  - (da  - Dc>*  - To  + fe:%  • 

The  three  boundary  conditions  for  the  third-order  system  (17)  are  seen 
to  be,  from  (14)  (13)  and  (16d) 

9(L)  - ♦(L)  - 0 . ♦ ’(<))«  0 . (18) 


With  ♦ ®ud  0 from  (17)  and  (18)  we  may  subsequently  determine  $ and 
v,  on  the  basis  of  (16a, b)  from  the  relations 


(19a) 

(19b) 


in  conjunction  with  the  boundary  conditions  $(L)  » v(L)  » 0. 

We  will  limit  ourselves  here  to  a determination  of  the  functions 
fandO.  The  value  0Q  of  0 at  the  loaded  end  z - 0,  when  written  in  the 
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form  0Q  - CqqQ0  + C0^.To,  then  gives  the  coordinate  xg  of  the  shear  center 
of  the  cross  section  In  the  form  -CRq/C0T  (3,4). 

We  find  from  Eqs.  (17a, b).  In  conjunction  with  the  second  and 
third  boundary  condition  in  (18),  as  expression  for  ip , with  D = DA+DB+D£, 


♦ - 

I <VDC>To  1° R 

Dc  2Da+D, 

b\ 

I ( i cosh  pz  \ 

1 4DAbC+DB ( DA+DC>  ' ^RM 

'4DaVDB(DA+DC) 

H1  cosh  pt/ 

where 

P2  - 

|4nADc  + 0b(Da  + Dc>I/DRD  • 

and  then,  writing  (17b)  with 

0 (L)  =0  in 

the  form. 

D0  - 

/'|(0A  - Dc)»  + T„  + 

<DC/DCQ>Qo,dz  . 

0 

V°C  [-  ( 

VDc  Dr 

Ml 

C 

o 

o 

0*= 

04 

- 1^ 

L 

*DADC+DB<DA+DC>’  ( ° ' 

da+dc  drm 

7 

VDc  f da-dc 

- 4bADc+DB(DA+Dc)  | c To 

(DR  2DA+DB  DC  ^ 1 sinh  pz-sinh  pL  (21) 

-\J5^ D Qoj  ITL'c'osra • 

Equation  (21)  as  it  stands  makes  it  possible  to  assess  the  effect 
of  finite  axial  length  L on  the  value  of  the  shear  center  coordinate. 

We  will  here  not  concern  ourselves  with  this  question  and  limit  attention 
to  cases  for  which  uL  is  sufficiently  large  to  allow  omission  of  the 
terms  with  p in  the  value  of  eQ.  We  then  get  from  Eq.  (21)  as  expres- 
sion for  the  location  of  the  shear  center 


x 


8 


^A~^C  2DA  DC 

VDC  ^RM  bA+DC  DCQ 


(22) 
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Having  (22)  we  may  determine  the  distance  of  the  shear  center  from 


the  elastic  centroid  of  the  cross  section,  as  follows.  We  have,  if  we 
designate  the  coordinate  of  this  centroid  by  x£ , E2A2(a  - xc)  = 

EjAjta  + xc)  and  therewith 


,E2A2'E1A1 


DR 

PRM 


(23) 


With  this  we  obtain  then  from  Eq.  (22) 


2D 


x - x 
s c 


A /DR  DC  \ 

W^rm  - JS^) 


CQ' 


(24) 


where,  from  Eqs . (8)  and  (9), 


DC 

^CQ 


,t2G2"tlGl 

t2(’2+tlGl 


(25) 


and 


1+f 


^tlGlt2G2 


(26) 


(tiG1+t2G2>ti2Gi2 


Equations  (23)  to  (25)  may  he  considered  to  be  the  principal  result 
of  the  present  considerations.  We  read  from  this  result  in  particular 
that,  Insofar  as  the  location  of  the  shear  center  relative  to  the  elastic 
centroidal  axis  of  the  beam  is  concerned,  this  location  depends  on  the 
distribution  over  the  cross  section  of  both  axial  normal  stress  stiffness 
and  transverse  shear  stress  stiffness.  At  the  same  time  this  location 
comes  out  to  be  independent  of  the  plate  twisting  stiffness  coefficient 
Dg  of  the  cross  section. 
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Two  limiting  special  cases  of  interest  are  given  by  (i)  the  case 
of  an  "open"  cross  section,  for  which  t^G^  - 0,  and  (ii)  the  case 
for  which  A^E^  =0,  as  shown  in  Figure  3.  We  have,  for  the  open  cross 
section  case  that  D^.  - 0 as  well  as  D^q  - 0,  in  such  a way  that  D^/D^q  * 
-a.  Equation  (22),  in  conjunction  with  (23)  then  reduces  to  the  form 


(O 


„ . E1A1-E2A2 

2a  + a 


(27) 


and  therewith  a<xa<3a,  depending  on  the  value  of  E^A^/E^j. 

When  AjE^  = 0 then  Dr  = 0 as  well  as  Dr^  = 0,  in  such  a way  that 
Dr/Drm  = -a.  We  now  have  from  (24)  and  (25) 


(ID  xs 


a 


—A—(l 

VDcV 


t2G2-t1G1\ 


(28) 


and  therewith -3a<xg<a,  depending  on  the  values  of  c/a,  t^G^,  t2G2  and 
Direct  asymptotic  determination  of  x . We  may  obtain  Eq.  (22) 
directly  by  considering  the  differential  equations  (17a, b)  for  the  limit- 
ing case  Dr  **  0 for  which  of  the  three  boundary  conditions  in  (18)  only  the 
condition  0(L)  - 0 remains  relevant.  In  order  to  obtain  the  correct 
result  for  this  limiting  case,  we  must  at  the  same  time  set  D^  - 0 
in  such  a way  that  remains  finite,  consistent  with  the 

relations  in  Eqs . (7).  W«  have  then  from  (17a) 


* 


VDC 


0 * 


(dr  dc  \ Q 


and  from  (17b) 


(29) 
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(30) 


(° 


(W 

a+db+dc  - -57^- 


)6’ 


t_  + q 


VDc  /dr 


4 % ^ V^RM 


•^K  • 


and  therewith 

*PADC+VDA+DC> 

da+dc 


(z  - L)  T 


PA~PC  PR 
*A+bC  °RM 


2Da  Dr 


(31) 


with  Eq.  (31)  once  again  implying  the  asymptotic  result  (22)  for  the 
quantity  xg. 

It  would  be  tempting  to  say  that  the  assumption  DR  » 0 means  that 
the  stress  measure  R becomes  a reactive  quantity,  that  is,  must  be 
determined  through  use  of  the  equilibrium  equation  R'  » T - T . 

w A 

However,  what  actually  happens  is  that  we  now  determine  R through  use 
of  the  constitutive  equation  (16d)  (with  DR  « 0 and  0^/D^  * 0)(  with 
subsequent  determination  of  a 7.- independent  expression  for  \|>  from 
the  equilibrium  equation  involving  R'  . 

A Remark  on  Generalizations . It  is  evident  that  an  analogous 
analysis,  involving  six  standard  force  and  moment  equilibrium  equa- 
tions, in  conjunction  with  a seventh  equilibrium  equation  of  the  form  (3d) 
can  readily  be  accomplished  for  a general  four-corner  cross  section  as 
depicted  in  Figure  A.  What  is  not  readily  seen  is  how  to  establish 
rationally  a seven  one-dimensional  equilibrium  equation  beam  theory  for 
a general  n-corner  cross  section  when  5£n  ( as  well  as  for  homogeneous 
cross  sections) , over  and  above  an  approximate  formulation  via  use 
of  the  Principle  of  Minimum  Potential  Energy,  in  conjunction  with  a suit- 
able four  parameter  family  of  expressions  (involving  St.  Venant's 
warping  function)  for  the  cross  sectional  distribution  of  axial  displace- 
ment, analogous  to  what  has  been  done  in  |A|. 
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One-dlmenslonal  beam  equations  are  derived  for  a rectangular-cross  section 
beam  consisting  of  shear  webs  and  corner  stringers  with  one  cross-sectional  axis 
of  symmetry.  It  Is  shown  that  a systematic  statement  of  this  problem 
Involves  one  more  equation  than  elementary  theory,  with  this  additional 
equation  involving  the  concept  of  cross-sectional  warping  and  bl-moment. 

As  an  application  of  the  general  results  of  the  paper,  new  conclusions  are 
deduced  concerning  the  location  of  the  center  of  shear  of  the  beam.  i—  
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